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Continuous melting of a driven two-dimensional flux lattice with strong pins 
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The phase diagram of a driven two-dimensional vortex lattice in the presence of dense quasi-point 
pins is investigated. The transition from the crystal to the liquid is found continuous at intermediate 
inductions. The correlations in the pseudo random force that allow for an uncomplete unbinding of 
the dislocations is proposed as a key mechanism to account for the continuous transition. 
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I. INTRODUCTION 

It has long been noticed that a driven elastic lattice 
driven at zero temperature may experience pinning as 
an effective shaking temperature, due to randomly in- 
duced displacements of the lattice nodes. This power- 
ful analogy allows for the prediction of some properties 
of the driven lattice from the common phase diagram 
of particles with a repulsive interaction. In particular, 
a dynamic melting transition is predicted and observed 
numerically jlj]. However, strongly disordered systems, as 
obtained in the presence of strong quasi-point pins, may 
move aside from this picture. Indeed, the thermal anal- 
ogy may break down for the moving crystal due to tem- 
poral correlations of the pseudo thermal Langevin force, 
a situation which is encountered in the case of hetero- 
geneous pinning involving plastic flow channels^- As 
a consequence, while this analogy accounts for the exis- 
tence of a first order-like dynamic transition, the driven 
phases may differ from their thermodynamic analogues. 
Several examples of such exotic phases have been given 
in refs. ^ A; 5, 6]. The anisotropy of the pinning poten- 
tial, once tilted by the external force, is essential to the 
formation of these phases0,l3l- Here, the dynamic tran- 
sition is examined in detail for a simple system of dense 
quasi-point pins with positional disorder. For part of the 
phase diagram, it is found that there is a continuous tran- 
sition between the crystal and the liquid, through what 
might be called a liquid crystal, whereas the transition 
is first order like for the rest of the phase diagram. In 
the case of the continuous transition, disclinations tend 
to form chains, which likely arise from the correlations of 
the pseudo random force which are specific to the driven 
lattice. 



II. EXPERIMENT 

The numerical sample used here is the one in ref. 
It is the one of two dimensional particles with a repul- 
sive interaction, interacting also with a random attractive 
potential. The sample mimics a two dimensional vortex 
lattice, or a three dimensional rigid vortex lattice, in the 
presence of strong pins, as can be created by heavy ions 
irradiation. Adopting the terminology of superconduc- 



tors, the vortex density is set by the magnetic induction, 
B, as n = Oq ^ = B/^q - $0 being the flux quantum car- 
ried by each vortex (2 10^^ Gcrri^). The repulsive force 
between vortices is taken as : 



UAr)^{Av/X)Ki (r/A) 



(1) 



where Ki is a Bessel function, behaving as In r 
at short distance and r^^/'^ exp(— r) at large distance. 
To keep computation tractable, the repulsive force is cut 
smoothly at a distance llao, which insures that each par- 
ticle interacts with many of its closest neighbors. 

The short range potential originates from strong pins 
randomly distributed in the sample, each creating the 
attractive force in the range rp: 

fp{r) = (2 Ap/rp) {r/rp), for r <rp; for r > rp (2) 

All pins arc identical and the randomness of the potential 
originates from the pins position only. 

In the rest, driving current densities are normal- 
ized to the single vortex critical current density, Jc = 
2 Ap/rp $0- The density of the pinning sites, relative 
to that of the vortices, B^/B, is constant and equal to 
12. The pinning potential range, relative to the vor- 
tex average separation, is also constant and equal to 
rp/flo = 5.5 10^^, as well as the reduced force mag- 
nitude, XAp /rpAy = 20. As a consequence, using oq 
as the length scale, the different numerical experiments 
made for different values of the induction B only differ 
by the reduced vortex interaction length, A/oq, where it 
was set A = 1400 A. As in the total force on each 
vortex, originating from its neighbors, a possible pinning 
site at the vortex location and the uniform external force 
is computed at each time step. Vortices which are not 
pinned are then moved on a time interval small enough so 
that their motion is small compared to all characteristic 
lengths. The boundary conditions are periodic along the 
driving force direction. A large area free from any pin- 
ning site is kept at the sample edges parallel to the vor- 
tex motion, where a perfect hexagonal lattice is obtained 
under the action of the external magnetic pressure*^. In 
this way, the measurements actually sample the driven 
phase embbeded in the crystal. Whereas such an inter- 
face may promote the formation of the ordered phase in 
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the case of a first order transition and for finite samples, 
in the case of a continuous transition, as will be discussed 
later, the interface probably induces an interfacial layer 
only. In the following, samples far enough from the inter- 
face are considered and their uniformity is an indication 
that finite size effects are not playing a major role when 
a continuous transition is observed. 

Experiments are carried out for different values of the 
induction and of the external force. After a stationary 
state is obtained (characterized by a steady average ve- 
locity), a snapshot of the moving lattice is recorded, on 
which a Delaunay triangulation is performed. Positive 
and negative disclinations (vortices with coordinance 5 
and 7), either free or forming dislocations by pairs 9J 
are counted. Samples typically enclose 7000 vortices and 
4 10"* pins. 



III. RESULTS AND DISCUSSION 

As shown in ref.|3, as the driving force decreases, 
the system evolves from a moving crystal to an amor- 
phous phase. Contrasting with the results in H, 
the high velocity phase does not show here smectic or- 
dering, as evidenced from the diffraction pattern : this 
comes from the small ratio rp/ao and from the fact 
that the tilted pinning potential shows here a moder- 
ate anisotropy on the scale of oq. There is no attrac- 
tive interaction between the vortices, which would allow 
for a transition between a liquid and a gas. However, 
considering the comparable densities of the crystal and 
the less ordered phase, as well as the strong interactions 
between the vortices in the amorphous phase, it must 
obviously be called a 'hquid phase'. As evidenced in 
FigH and El and the inspection of the average hexatic 
parameter,! < ^-g > | = | < 1/ca E/3=i,c„ * >" I 
where is Cq. the coordination number for vortex a and 
Sa,i3 is the angle of the bond between neighboring vor- 
tices a and /3, some residual orientational correlation is 
retained for low j (| < vj/g > | ~ 0.1), which justifies to 
call the low j phase an 'hexatic liquid' Q. 

I now examine in more detail the transition between 
the crystal and the liquid. For all systems, the concen- 
tration of defects exhibits a clear onset upon decreas- 
ing the driving force, similar to the one reported in 0. 
However, depending on the magnetic pressure, a dis- 
continuous or gradual rise of this concentration is ob- 
served. This may be seen in Figs. H and El obtained 
for two different magnetic inductions, which clearly ex- 
hibit respectively a gradual and a step increase of the 
number of defects. In order to quantify this observation, 
the defects concentration was fitted with an exponential, 
cx 1 - exp[(jo - j)/6] {j < jo), yielding the onset, jo, 
and a width for the transition to the liquid phase, S. A 
phase diagram similar to the temperature-density repre- 
sentation for the thermodynamics may be obtained, using 
the theory for the equivalent 'shaking temperature' [ij. 
It should be stressed that this representation is qualita- 



tive only, considering the reservations made in ref . jlj, |3| 
(mainly, the perturbative approach from the uniform ve- 
locity which rules out plasticity, and the observation that 
the effective temperature differs for the fluidlike motion 
and the coherent one). Also, the equivalent temperature 
in 1], must be modified to account for the proximity of 
the flux flow to the flux creep crossover. Imposing for 
the equivalent temperature to be proportional to the po- 
tential well depth when io — > 1 and B— >0,Tcx(l — j), 
and using T oc from |l|, a phase diagram is obtained 
in the B vs (1 ~ j)/j representation. The onset for the 
defects concentration (jo), as well as the location where 
it saturates {jo — S), are plotted in this way in FigO 
Clearly, there is a range of magnetic induction for which 
a regime, intermediate between the moving crystal and 
the hexatic liquid, can be found. The existence of such a 
regime was already pointed out in 8] . 

In order to characterize the continuous transition, let 
us examine some autocorrelation functions which are 
classical tools for the study of solids and liquids. The 
average hexatic order parameter does not provide an ac- 
curate characterization of the intermediate regime: as 
may be seen in Fig. ^ following a sharp drop at j — jo, 
there is no significative change at lower j where the 
density of defects however still exhibits significant vari- 
ations. The spatial correlations of the hexatic param- 
eter carry more useful information 'o'j. The correlator 

< ^'g(O) ^'g(r) >r for the data in Figgis shown in 
FigEl Besides the existence of a non zero background 
related to the non zero averaged value | < ^fg > |, it 
reveals some additional short range correlations of the 
orientational order, which extends to a larger range as 
the system is closer to jo- The oscillations for small r re- 
flect the existence of a crystalline order within this range: 
they are associated with the fluctuations of the density 
autocorrelation function which come with the transla- 
tional symmetry breaking of the crystal order. This is 
confirmed by the examination of the growth of the dis- 
placement field (actually a positional correlation func- 
tion): < u^(r = n Co) >=< J2j u'^in [rj — r^]) >i where 
j denotes one of the nearest neighbors of vortex i in a 
Delaunay triangulation and u is the displacement field 
from the periodic arrangement. As may be seen in Fig. 
there is an exponential decay of the positional cor- 
relations, with a diverging correlation length, £,{j), as 
one approaches jo from below (i.e. from larger 'tempera- 
tures'). Following ref. j^, one may then call this regime 
an 'hexatic liquid crystal'. It is possible to track the po- 
sitional correlation length in FigElas one approaches jo- 
The result is diplayed in FigO showing a divergence as 
C = ^0 (1 ~ j/jo)^^, with the bare correlation length 
^0 — 0.15 flQ. Recalling the Lindeman melting criterion, 

< >= c| Oq and the exponential increase of the dis- 
placement field, < >=< u^(oo) > (1 — exp(— r/^)), 
one may write an equivalent melting criterion for the 
present case as ao/^ = ln(l — c| Oq/ < u^(oo) >). The 
result obtained using cl ~ 0.2 and < u'^{oo) >~ 0.14, 
oo/C ~ 0.3, is displayed in Fig^andEl Although this 
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quantitative result should be considered with caution, 
due to the uncertainty on the effective Lindeman num- 
ber, this confirms that the solid has not melted in the 
conventional way below the threshold value jo- 

The observation that the transition is continuous at 
intermediate induction and involves a proliferation of de- 
fects appeals for a comparison with the KTHNY exten- 
sion of the Kosterliz-Thouless theory ^| . The theory 
of dislocation mediated melting of two-dimensional solids 
accounts for a continuous transition, involving first un- 
binding of the dislocations leading to the hexatic liquid, 
and then unbinding of the disclinations leading to the reg- 
ular liquid. Here, dislocations do not first dissociate at ja 
to form an homogeneous 'plasma'. Rather, they tend to 
form chains of alternating positive (five-coordinated) and 
negative (seven-coordinated) disclinations which prolifer- 
ate in the liquid phase. As a result, unbounded disloca- 
tions and disclinations remain marginal (Fig. EJ. Corre- 
lations between dislocations were also reported in ref. 4] 
where free dislocations, although not bound in chains, 
formed quenched patterns moving with the average flux 
flow. In order to explain the formation of these chains, 
the examination of the early creation of defects in a 
driven crystal may be useful. Snapshots of the earlier 
defects detected in a sample driven in the intermediate 
region in Fig^are displayed in FiglHl After a dislocation 
pair with opposite Burger vectors has been created by 



the pinning of one vortex (a), it is seen that the dislo- 
cations quickly arrange to form rings of diameter ~ 2 ag 
(c) and then larger loops (d) . Remarkably, the composite 
defects reflect the external force anisotropy as soon as the 
dislocations dissociate (b): this results from the plastic 
mechanism at work to create these defects. This is also 
a direct evidence that the correlations in the pseudo ran- 
dom force cannot be neglected in their formation. The 
relation between these initial stages and the formation of 
chains is not completely clear. A possible mechanism is 
the stretching of elementary loops as in FiglHli, as the vor- 
tices making disclinations appear to become more easily 
pinned than the regular ones. This would make the long- 
range 'random force' correlations a key ingredient in the 
chain formation again. Equivalent rates for the growth 
and the annihilation of the chains would then account for 
the existence of a stationary regime intermediate between 
the crystal and the liquid. 

In conclusion, it is found that a vortex lattice driven 
on dense quasi-point pins shows a continuous transition 
between the crystal and the liquid at intermediate in- 
duction, while first order otherwise. The binding of the 
disclinations in chains is proposed as a key mechanism to 
account for the existence of the continuous transition. 

Simulations have been performed on the cluster of the 
Centre de Ressources Informatiques de I'Universite Paris- 
Sud (CRI). 
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FIG. 1: B = 3950 Oe. From top to bottom : Hexatic parame- 
ter, concentration of defects (sites with coordination number 
not equal to 6), concentration of free disclinations (defects 
bound to sites with coordination number 6 only). The line is 
the fit described in the text; the full line arrow indicates the 
onset for the defects creation, as obtained from this fit; the 
dotted one is the melting point as obtained in Fig|S| 




FIG. 2: B ^ 10* Oe. 
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FIG. 3: Dynamic phase diagram of the driven lattice. Circles: 
j = jo, triangles : j = jo - 5. 
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FIG. 4: The correlation function for the hexatic parameter. 
{B = 3950 Oe) 
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FIG. 5: The correlation function for the displacement field. 
{B = 3950 Oe) 
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FIG. 6: B = 3950 Oe. Correlation length for the positional 
order, as obtained from the data in Fig|3 The full line is a 
linear fit. The dotted line represents < u'^ >= ao with 
cl ~ 0.2. The melting point, as indicated by the dotted 
arrow, is shown in Fig also. 




FIG. 7: Positive - five-coordinated (white) and negative - 
seven-coordinated (black) disclinations in a sample driven 
along the vertical axis . {B — 3950 Oe, j = 0.375) 



FIG. 8; Early steps for the creation of defects {B = 3950 
Oe, j = 0.375. Nodes of the Delaunay triangulation hnes 
are vortices. Filled and opened circles are disclinations as 
in Fig|7| gray circles indicate vortices which are located in a 
potential minimum). Configurations a) and b) are obtained 
consecutively as the result of the trapping of one vortex. Con- 
figuration c) and d) are later steps. Vortices are driven to the 
top of the figure. 



